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ABSTRACT 

We present 838 afe-type RR Lyrae stars from the Lowell Observatory Near Earth Objects Survey 
Phase I (LONEOS-I). These objects cover 1430 deg 2 and span distances ranging from 3-30kpc from 
the Galactic Center. Object selection is based on phased, photometric data with 28-50 epochs. We use 
this large sample to explore the bulk properties of the stellar halo, including the spatial distribution. 
The period-amplitude distribution of this sample shows that the majority of these RR Lyrae stars 
resemble Oosterhoff type I, but there is a significant fraction (26%) which have longer periods and 
appear to be Oosterhoff type II. We find that the radial distributions of these two populations have 
significantly different profiles (p 0o i ~ R~ 2 - 26±007 and p 0o ii ~ i?" 2 - 88±011 ). This suggests that the 
stellar halo was formed by at least two distinct accretion processes and supports dual-halo models. 
Subject headings: galaxies: structure — galaxies: halos — galaxies: individual (MWG) 



1. INTRODUCTION 

A major goal of modern astrophysics is to understand 
the physical processes by which galaxies form, and to 
place galaxy formation in the context of models for the 
evolution of large-scale structure in the Universe. The 
advent of wide-field CCD cameras has now made it pos- 
sible to survey large samples of Galactic stellar popula- 
tions, and therefore to study the structure of our Galaxy 
as a whole. While keeping in mind that the Galaxy may 
not always be representative of (spiral) galaxies in gen- 



eral, the study of our Galaxy is now a fruitful avenue for 
the detailed exploration of galactic structure, providing 
the empirical basis for differentiating between formation 
models. 

RR Lyrae stars are horizontal branch stars which lie 
in the instability strip. They represent a late stage of 
stellar evolution, and as such form an old stellar popula- 
tion that can be used to study the early stages of galaxy 
formation. RR Lyrae stars have a long history as tracers 
of Galactic structure ( Kinman et al.lll96ll lOort &: PlautJ 
[19751 IHawkinsI fl98l ISahal I1985D Thev are also rela- 
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tively easy to identify from their colors and variability, 
enabl ing the format i on of s amples with little contamina- 
tion. ISuntzeff et alj (|1991h estimate the number of field 
RR Lyrae in the Galactic halo to be about 85,000, but 
less than 8,000 have been dis covered, inc luding those in 
the bulge and thick disk (e.g., Heck 1988). Further, their 
absolute magnitudes (for RR Lyrae afe-types) have been 
measured to bett er than 0.1 magnitu des using robust sta- 
tistical para llax (lLavden et al.lll996l) and direct parallax 
techniques (|Benedict et al.ll2002l) . However, RR Lyrae 
stars do have a disadvantage. The relative abundance 
of RR Lyrae stars depends on the horizontal branch 
morphology of the RR Lyrae progenitor system. Thus, 
they may not reveal the entire history of Galactic for- 
mation. Overall, RR Lyrae stars are an incisive probe 
of the halo field population and may provide useful in- 
sights on Galactic formation. In recent years, several 
sky survey have discovered large numbers of Galactic 
RR Lyrae oyer la rge fractio ns of the sky e.g, SDSS 
(llvezic et alJ[2000h. QUEST (I Vivas fc Zinni 120061 and 
NS VS (jKinemuchi et all [20M) . The SDSS and QUEST 
surveys have shown that RR Lyrae can be used to find 
and trace the remants of accreted satellite galaxies in the 
Galaxy. 

RR Lyrae also have been well cataloged and stu died 
in the globular clusters fe.g.. IcTement et al.ll200lT ). In 
particular, globular clusters fall into two broad groups 
based on the mean periods of their RRab stars and on 
their RRc fraction, with Oosterhoff I (Ool) clusters hav- 
ing shorter mean periods and lower RRc content com - 
pared to Oosterhoff II (OoII) clusters (|Oosterhofi1ll939D . 
In addition, these two classes populate distinct regions 
in period-amplitude space. There are certain stellar sys- 
tems (most Local Group dwarf spheroidals) which can- 
not be classified as Ool nor OoII, but rather exhibit 
period-amplitude distributions which fall between Ool 
and OoII. These are referred to as intermediate. We 
can examine Galactic stellar halo RR Lyrae in a simi- 
lar manner to deter mine their Oosterhoff classification. 
ISuntzeff et alj (1991) found Ool and OoII components in 
field RR Lyrae stars, but no intermediate Oo component 
(i.e., the Oosterhoff "gap"). The origin of this Oosterhoff 
effect has been studied intensively (e.g., iSandage et al.l 
[l98ltlLee et "all [19951 ILee fc Carnevlfl99l . The Ooster- 
hoff classification may either indicate an age difference 
or a metallicity difference. Regardless, the Oosterhoff 
classification of the RR Lyrae in the Galactic stellar halo 
will provide insights into it formation history. 

In this paper, we present new results from the LO- 
NEOS - I RR Lyrae survey. Using 838 RR Lyrae stars 
discovered from the LONEOS - I survey, we will investi- 
gate their distribution in the stellar halo out to distances 
of 30 kpc, as a function of Galactic position. In Section 2, 
we introduce the LONEOS - 1 variability survey and dis- 
cuss the photometric pipeline. In Section 3, we present 
the LONEOS - I RR Lyrae sample and discuss the se- 
lection process and detection efficiency. In Section 4, 
we discuss the properties of the LONEOS - I RR Lyrae 
stars. Finally, we discuss the implications for Galactic 
formation and look to the future. 

2. LONEOS - I VARIABILITY SURVEY 

The Lowe ll Observatory Nea r Earth Object Survey 
(LONEOS) (jBowell et al.l)l995f) has as its primary goal 
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Fig. 1. — The black regions indicate subfields of LONEOS - I 
survey data with at least 28 good quality images. These fields are 
used for our RR Lyrae star analysis in this paper. This area covers 
1430 deg 2 on the sky. The equatorial coordinates of the subfields 
are displayed using a Mollweide projection. 

the detection of potentially hazardous asteroids. Situ- 
ated on the Anderson Mesa outside Flagstaff, Arizona, 
LONEOS is carried out with a 0.6m Schmidt telescope. 
During the first two years of operation, imaging was done 
with a camera construct ed at the University of Washing- 
ton (|Diercks et al.l ll995). with a single 2k x 4k unfiltered 
CCD (5 sq. deg. FOV with 2.67"/pixels) (LONEOS - 
I ). Under optimal conditions, LONEOS can image up 
to 1000 square degrees per night, multiple times. The 
LONEOS - I survey imaged fields three times a night 
with 60 second exposrure times. Fields are re-observed 
at least one other time about one lunation later. LO- 
NEOS - I camera uses no bandpass filter, but the data 
can be transformed to a standard photometric system 
using external datasets. The LONEOS - I camera can 
reach a depth of R ~ 18.5. 

Besides searching for moving objects, this dataset can 
be used to search for temporal variability, specifically of 
stellar objects. Data from LONEOS have been used to 
construct a stellar temporal variability database, from 
which we extracted candidate RR Lyrae stars. The ca- 
dence of LONEOS observations makes them well suited 
to detect short-term variability (e.g., RR Lyrae, binaries, 
CVs), and long-term variability (e.g., QSOs). 

The LONEOS - 1 photometric variability database dis- 
cussed in this paper consists of photometric data ob- 
tained between 1998 - 2000. It covers a large fraction 
of the sky, much of it imaged in multiple epochs (see 
Figure [1]). The LONEOS - I dataset that we will ana- 
lyze in this paper consists of 1430 deg 2 with at least 28 
epochs. 

2.1. Image Reduction Pipeline 

The LONEOS - I image reduction pipeline consists of 
several steps. Basic CCD processing (bias, flat-field cor- 
rection) is followed by object detection using the photom- 
etry program DoPhot (|Schechter et al.lll993D . Bookkeep- 
ing is accomplished by automatic merging of the DoPhot 
output file with the FITS header of each parent image to 
produce a self-describing and self-contained output file. 
These output files are entered into a database (broken 
down by region on the sky) where the primary index- 
ing is based on image position (RA, DEC). The relative 
and absolute photometry is performed at the database 
level. We divide each field into subfields, and relative 
photometry is performed on these subfields. We tie the 
LONEOS - I absolute photometry to an external cata- 
log. Details are presented in later sections. Access tools 
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are used to search the database and apply variability 
detection algorithms. We have developed a set of ro- 
bust statistics which are used to discriminate between 
long and short term variables in the timeseries data. It 
should be noted that we do not attempt to perform ac- 
curate photometry of extended objects. The database 
should be strictly used for stell ar ob j ects. For details of 
the data reduction pipeline see iRestJ (|2002D . 

2.1.1. Photometry 

The photometric calibration can be separated into two 
independent parts: relative photometry and absolute 
photometry. For identification and classification of vari- 
able objects, relative photometry suffices. Since LO- 
NEOS uses an unfiltered system, we must rely on exter- 
nal datasets to tie the relative magnitudes to a standard 
photometric system. We outline the procedure in this 
section. 

For an image taken at a given time t and airmass a(t) 
through a filter X, the apparent magnitude m(t, s) of a 
star s with color X — Y (e.g. X — Y = B — V) can be 
written as 

m(t, s)=mi nstr (t, s) + A x (t) + c\ * a(t) 

+c 2 * (X — Y) * a(t) + c 3 * {X - Y) + m (l) 
Here mi ns t r (t, s) is the instrumental magnitude with an 
error of C7i ns t r {t, s), Ax(t) is the extinction due to clouds 
(assuming 'grey' extinction), c\ and C2 are the grey and 
color-dependent airmass extinction coefficients, respec- 
tively. For our purposes, since we are interested in the 
photometry of only one type of star, we can set (X—Y) to 
a constant, and include C2 in c\ and C3 in rriQ. LONEOS 
takes images of fields with generally low airmass (76% of 
images have airmass less than 1.5 with a mean airmass 
of 1.3). Therefore, the second order term depending on 
color and airmass is very small anyway. We group all 
image-dependent terms into m ca i(t) — Ax(t) + c\ * a(t) 
and 

m(t, s) = m mstr (t, s) + m ca i(t) + m . (2) 
A wide field-of-view or strong PSF distortions of the op- 
tical system can cause systematic variations in the pho- 
tometry across the image. We divide each LONEOS field 
(i.e, 2k x 4k) into 50 square subfields, each 0.304 degrees 
on a side. This subfield size is large enough to include 
a sufficiently large number of stars to perform the pho- 
tometry, but small enough to avoid systematic problems 
(e.g., Flat fielding errors, clouds, etc) across the wide 
FOV of the camera. 

2.1.2. Relative Photometry 

We define the relative magnitude m re i(t,s) and the 
time-averaged relative magnitude rh re i(s) of the star s 
as: 



star measurements to the time-averaged star magnitude 
is minimized. The function that is minimized is: 



m re i(t, s) =m inst r(t, s) +m ca i(t), 
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We take the following approach to perform relative pho- 
tometric calibration: calculate the time-averaged mag- 
nitude for each star, then adjust m ca i(t) in each sub- 
field, so that on the average the deviation from single 
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where m ca i(t) are the free parameters. Note that without 
the second term, the x 2 is degenerate to the addition of a 
constant to the m ca i(t) vector. Therefore we add the sec- 
ond term which is quadratic in the degenerate direction. 
Hence our solution has ^2 t m ca i{t) — 0. We also add 
0.02 magnitudes in quadrature to <Xj„ sfrr (t, s) to take into 
account systematic error. This is intended to compen- 
sate for the fact that DoPhot is known to underestimate 
photometric errors at the bright end, where the error 
is dominated by systematics rather than Poisson noise. 
These systematics include flat fielding errors, PSF varia- 
tions, and the limited accuracy of the PSF model. Stars 
for which the reduced x 2 °f their lightcurve exceeds 2.5 
are flagged as potentially variable and excluded from the 
set of stars used to fit m ca i(t), with 

xl = l/(N t - l)J2(rh re i(s) - m rel (t, s)) 2 /a 2 {t, s).(8) 
t 

Here N t is the number of images. Similarly, we exclude 
images for which their reduced xf exceeds 1.5, with 

X 2 - l/(N s - 1) Y,{fh rel (s) - m rel (t, s)) 2 /a 2 (t, s).(9) 

s 

In addition, we exclude images for which more than 20% 
of the star measurements m re i (t, s) deviate by more than 
3 a from their time-averaged magnitude m re i{s). This 
process is repeated iteratively until the m ca i(t) vector 
values for each subfield do not change significantly. 

2.1.3. Absolute Photometry: Tying LONEOS Photometry to 
External Datasets 

Since LONEOS uses an unfiltered camera, we 
must use external datasets to tie LONEOS photom- 
etry to a standard photometric system. There exist 
datasets of calibrated photometry over large por- 

tions of the sky. We use the Guide Star Catalog 2.2 

(ISpace Telescope Science Institute k, Osservatorio Astronomico di r . 

2001) (hitherto, referred to as GSC). GSC is complete 

to 18.5 photographic F-band, and Tycho 2 photometry 

is used for bright sources. Thus, the GSC magnitude 

range matches well to LONEOS - I . The procedure 

to tie LONEOS relative photometry to GSC is similar 

to the relative photometry algorithm. We restrict the 

range of color (—1.1 < (Mgsc.j — Mgscf) < 1.5) for 

the GSC stars to be similar to colors of RR Lyrae, which 

are the focus of this work. The spectral response of the 

LONEOS - I camera system matches closest to the GSC 

F passband. Initially, LONEOS - I photometry is tied 

to GSC F, which we refer to as Mloneos_f with error 

&LONEOS-F- Combining Equations O [3] and [H we define 

the time-averaged apparent magnitude MloneoS-f(s) 

of the star s in each subfield t as: 



M L oneoS-f(s) : 
a(s)- 



- ~2 



ei{s) +m QiGSC _ F , (10) 
iW+Ccscf ( n ) 
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Fig. 2. — The residual between LONEOS - I magnitudes (tied 
to GSC2 F) and V S dss as a function of SDSS color, g - r. The 
two dashed vertical lines indicate the color range typical for RR 
Lyrae. In this range, we see that there is an offset, M^qneoS-F — 
V S DSS = -0.32, with an RMS of 0.13. 



Fig. 3. — The residual between LONEOS magnitudes (tied to 
GSC2) and SDSS as a function of RA for equatorial region. One 
can see the systematic difference between different POSS plates. 



We determine the free parameter mo^GSC-F and its error 



"mo,GSC-F 



for a given subfield by minimizing 



Right Ascension. One can clearly see systematic errors 
as a function of RA. These are on the scale of Schmidt 
plate size (6 degrees). However, these large errors occur 
over only a small fraction of the sky. 0.98% of the stars 

E, ., _ 2 have residuals greater than 3 a from the mean compared 

(MascM*) ~ MloneoS-fW) /fec.f + ^1% 0.3% for a Gaussian distribution. Since the survey 
s area which we will analyze is large, these errors do not 

affect the statistical properties of our RR Lyrae sample. 



where the sum is taken over the stars in a given image. 
Once mo is dctercmincd for each subfield, the LONEOS 
magnitude tied to GSC F for a given star s, in that sub- 
field is: 

MloneoS-F (t, s) (t, s) + m cal (t) + mQ.GscQS) 

However, the absolute magnitude of RR Lyrae have 
been measured in V. We use photometry from SDSS 
Data Release 3 (DR3) (jAbazaiian et al.ll2005l ) to convert 
MloneoS-F to a U-based system, in this case a syn- 
thetic SP SS-based magnitu des (Vsdss — r + 0A4(g — 
r)-0.02 (|Ivezic et al.l l2005) ). All magnitudes have been 
correc ted for Galactic extinction using Schleg el et al.1 
(1998). The residual of M LO neoS-F and Vsdss as a 
function SDSS color, g — r, has a clear trend (Figure [5]). 
Thus, we must restrict the range of SDSS color used to an 
empirical offset between Mloneos^f and Vsd££- Stars 
in the DR3 with colors of RR Lyrae (jlvezic 
are matched to LONEOS - I stars (0.99 <u- g < 1.28, 
-0.11 < g - r < 0.31, -0.13 < r - i < 0.2, -0.19 < 
i-z < 0.23). A total of 40,180 SDSS DR3 stellar ob- 
jects in this color range were matched to LONEOS - I 
objects. Using LONEOS - I timeseries, we reject stars 
with signs of variability. We see that there is an offset, 
MloneoS-F — Vsdss = —0.32, with a scatter of 0.13. 
All subsequent LONEOS - I photometry in this work is 
tied to the Vsdss system via this offset. The final trans- 
formation is: 

Mloneos v — [Mloneos f— Agscf xE(B— V)]+0.32), 

where Aqscf = 2.65 is the extinction coefficient and 
E(B-V) is the reddening. This transformation is only 
valid for stars with colors similar to those of RR Lyrae 
stars. 

The distribution of the residuals is approximately 
Gaussian in the core, but there are long tails in the dis- 
tribution, most likely due to plate-to-plate systematic er- 
rors in the GSC photographic sky survey. Figure [3] shows 
the residual of Mlqneos-F and Vsdss as a function of 



3. LONEOS - I AB-TYPE RR LYRAE STAR CATALOG 

In this section, we will describe the selection of RR 
Lyrae using LONEOS - I data, and present the LO- 
NEOS - I RR Lyrae star catalog. This currently rep- 
resents the largest single sample of halo RR Lyrae stars. 
In addition, we discuss the detection efficiency of our se- 
lection process and the completeness of our sample. 

3.1. Selection of LONEOS - I RR Lyrae Stars 

The selection of the RR Lyrae was solely based on LO- 
NEOS - I temporal photometric data. No color informa- 
tion was used. We restrict our search to fields with more 
than 28 good quality images. This amounts to 1430 deg 2 
of coverage. The median number of images is 31. We 
impose a set of variability criteria which first identifies 
short-term variable stars, and then analyzes the morphol- 
ogy of the phased light curves. We limit our search to 
aft-type RR Lyrae, since it is difficult to detect efficiently 
the c-type RR Lyrae given our temporal sampling. In 
order to accurately determine our detection efficiency, 
we did not directly rely on visual inspection to produce 
our sample, but used a series of cuts which aim to maxi- 
mize the detection of RR Lyrae in an objective manner. 
Below, we outline our selection criteria (see Figure 2]). 

To reject non-variables and/or low amplitude variables, 
we made two variability cuts on the timeseries: 

• We first looked at the spread of the measured 
magnitudes about the average magnitude, nor- 
malized by the photometric error, which is a re- 
duced x 2 - We define this statistic as: x 2 = 
Y^riLi ( m '~/") , where N is the number is ob- 
servations, m, is the magnitude of the i th obser- 
vation, m is the mean magnitude, and cti is the 
photometric error of the i th observation, x 2 gives 
a measure of the variability over the entire time- 
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Fig. 4. — LONEOS - 1 RR Lyrae flow chart indicate the selection 
process and criteria used in the analysis of this paper. 

series. We rejected stars with \ 2 < 3-5. This pre- 
dominantly rejected small amplitude variables. 

• We divided the timeseries into two segments, each 
containing half the data points, and calculated % 2 
for both halves. We required x 2 > 2.5 for both 
halves of the timeseries. This rejects objects that 
are aperiodic or that have inflated \ 2 from occa- 
sional outliers. 

A total of 5889 objects satisfied these cuts. In order to 
identify RR Lyrae stars in this sample, we search for 
a periodic signal in these timeseries. In addition, we 
must classify Period-folded light curves. However, the 
LONEOS - I sampling is not sufficient to use sophisti- 
cated classification methods, such as Fourier decomposi- 
tion or machine learning methods. We use template light 
curves to distinguish RR Lyrae from other common light 
curves, such as sinusoidal variab les and eclipsin g binaries. 
We use the nine templates from lLavd cn (1998), based on 
high quality light curves of bright RR Lyrae stars which 
have been averaged. Six templates are ab-type RR Lyrae 
stars. In addition, there are single templates for c-type 
RR Lyrae, for a sinusoid, and for a W UMa-type con- 
tact binary (Figured]). There are a sizable number of 
contact binaries in the LONEOS - I database. The fun- 
damental mode RR Lyrae afe-type can be broken down 
into two different categories, type a and b (i.e., "Bailey" 
type). Bailey type a objects have flatter bottoms, and 
the Bailey type b have a sawtooth shape. The first six 
templates represent typical RR Lyrae light curve shapes 
for type a and b with three different rise times (i.e., the 
phase interval between maximum an d minimum ligh t). 

We used a supersmoother routine (Rc imannlll994f ) to 
obtain a set of four best periods. The supersmoother is 
a variable span smoother, where the best span is chosen 
using a cross-validation te chnique. It was de veloped for 
the MACHO experiment (jAlcock et al.lll995h . We have 
tested other st andard routines, su ch as Phase Dispersion 
Minimization ( Stcllingwcrf 1978) and Lomb-Scargle pe- 
riodogram (|Scargld 11981 : these alternative algorithms 



Fig. 5. — The best four periods from the supersmoother are fit to 
these nine templates. The mean magnitude, phase and amplitude 
are free param eters in the fit. These nine templates are taken from 
ILavdenl 1119981 ). The first six are afe-type RR Lyrae which were 
obtained from the averaging of high-quality light curves. The last 
three include a c-type RR Lyrae, a sine wave, and W UMa contact 
binary. 



provide comparable results to the supersmoother. The 
supersmoother makes no assumptions regarding the mor- 
phology of the light curve. As a consequence, the super- 
smoother's best period may not represent a physically 
realistic light curve. In addition, we have found that 
light curves with outlying data points adversely affect 
the supersmoother's figure of merit for smoothness. 

We limited the period search to range between 0.25 
and 0.9 days. Our results were not affected if the search 
is extended to 2.0 days. In order to break degeneracies 
between the four best periods returned, we used a non- 
linear \ 2 minimization method to fit the phased data to 
our set of template light curves (Figure [5]) . The magni- 
tude offset, phase and amplitude were free parameters in 
the template fitting. We required that the fit amplitudes 
be greater than 0.3 magnitudes. Each of the four periods 
was fit separately. The best period was selected by the 
quality of the template fit. 

Objects were tentatively classified as afe-type RR Lyrae 
if convergence of the non-linear minimization was satis- 
fied for the fits to one of the six RRafr templates. We 
defer any template light curve goodness-of-fit cuts un- 
til later. Objects with best fits to the three non-RRa& 
template light curves were discarded. A total of 1238 ob- 
jects were classified as a&-type RR Lyrae. These objects 
were visually inspected and classified into three groups: 
"good", "fair", and "bad". The criteria for the classifi- 
cations included: goodness-of-fit to template light curve, 
and uniformity of data across phase. It is fairly easy to 
classify objects in the "good" and "bad" categories. The 
"fair" group includes objects where there was some am- 
biguity in the classification. These classifications were 
used as a training sample for further selection criteria. 
The aim was to determine a totally automated set of se- 
lection criteria to maximize the selection of real RRa& 
while minimizing contamination. This will facilitate the 
calculation of our detection efficiency. We developed the 
following criteria for quality of the light curves. First, 
we required that the fitted periods and amplitudes re- 
side in the region of period-amplitude space where RRa& 
are known to lie. This region was defined to contain the 
bulk of our "good" sample, and was large enough to en- 
compass the range of metallicities and Oosterhoff groups 
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associated with halo RR Lyrae stars. We are guided by 
the MACHO RR Ly rae period-amplitud e distribution for 
the Galactic Bulge (jAlcock et al.lll998t ). which contains 
a population of RR Lyrae stars with a wide range of 
metallicities. The region in period-amplitude space was 
defined as: 

A < 1.5 magnitudes, (15) 

A > 0.3 magnitudes, (16) 

A > —5 log(P) — 1 magnitudes, (17) 

A < -5 log(P) + 0.3 magnitudes, (18) 

where periods are measured in days. Using the well- 
sampled light curves from Nort h ern S ky Variability Sur- 
vey (NSVS), iKinemuchi et alJ (|2Q06f > found a popula- 
tion of a6-type RR Lyrae with shorter periods than 
the Ool group. Using light curve Fo urier decomposi- 
tion techniques. IKinemuchi et al.l (|2006l ) determined that 
this population was more metal-rich than the rest of 
the ROTSE-I sample. This population was identified 
as a metal-rich thick disk RR Lyrae star population. 
Our period-amplitude cut largely excludes this metal- 
rich thick disk RR Lyrae star population. A metal-weak 
thick disk RR Lyrae star population also likely exists, 
but this population smears with the Ool in the period- 
amplitude diagram. Thus, we are confident that out our 
RR Lyrae star sample will consist predominantly of halo 
objects. Of course, kinematic information is required to 
be certain of population membership. Finally, an OoII 
component may extend to smaller amplitudes than an 
Ool component. The amplitude cut should optimally 
be extended to smaller amplitudes at longer periods, but 
our detection efficiency below 0.3 magnitudes is relatively 
low at the faint end. While this cut may preferentially re- 
move OoII RR Lyrae stars, the LONEOS - 1 dataset can- 
not efficiently probe this region period-amplitude space. 

Using our visual classifications (i.e., "good", "fair", 
and "bad"), we determined the thresholds for two addi- 
tional cuts: uniformity and smoothness. These cuts were 
performed simultaneously, and their thresholds were set 
to minimize the number of "bad" objects and maximize 
the number of "good" objects. First, we required that 
the light curves have a degree of uniformity. The phased 
data points should be distributed from zero to one as ran- 
domly and uniformly as possible. In order to measure 
this, we divided the phased light curve into 10 equally 
spaced bins, and binned the data points. We required 
that at least 7 bins have at least one measurement. Thus, 
at most 30% of the phase is allowed to contain no data 
points. We referred to this as the "N gap" test. Secondly, 
we examined the goodness of fit of the RR Lyrae star 
template relative to the fit to a constant magnitude. We 
define this statistic as: A X 2 = {xlonstant ~ Xrr)/Xrr, 
where we use a reduced \ 2 f° r au the fits. This statis- 
tic measures the scatter about the six RRafr templates, 
and thus smoothness of the period-folded data. For this 
sample, we used Ax 2 > 5. 

Together, these three cuts reject 8%, 60%, 97% from 
the "good", "fair", and "bad" samples, respectively. Fig- 
ures [5] and [7] show the results of the final three cuts. The 
sample presented in this paper are restricted to mean 
magnitudes ranging from 13.0 to 17.5. We will discuss 
the completeness of this sample in the following sections. 
We also removed known RR Lyrae that are members of 
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Fig. 6. — These three period-amplitude distributions show the 
objects which were successfully fit to an RR Lyrae template light 
curve. The plots show the amplitude (in LONEOS - I instrumental 
magnitudes) versus the logarithm of the period measured in days. 
They have been manually classified into three categories, "good", 
"fair", "bad". The "good" and "fair" distributions are consistent 
with being RR Lyrae with some scatter, while the "bad" category 
is randomly distributed. The dashed lines indicates the region 
of period-amplitude space where RR Lyrae stars are expected to 
lie. This region is large enough to encompass a wide range of 
metallicities and Oosterhoff types. 
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Fig. 7. — We show the resulting objects after we have applied the 
period-amplitude cut and two additional light curve quality cuts 
described in the text. The lower-right panel shows the final RR 
Lyrae sample. 

Period=0.47 Amp=1.0 Period=0.63 Amp=0.8 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

phase phase 

Fig. 8. — Two examples of phased RR Lyrae light curves. The 
lightcurves are shown over two full phases. 

globular clusters and RR Lyrae within the tidal radius 
of the globular clusters in our survey volume, a total of 
19 objects. Our final sample consists of 838 aft-type RR 
Lyrae, which are solely based on temporal photometric 
data. Table [T] lists the properties of this sample. Two 
sample light curves are shown in Figure [8] To date, this 
represents the largest single sample of period-folded RR 
Lyrae stars in the stellar halo. 

As a simple test of the contamination of our RR Lyrae 
sample, we explored the colors of a subset of these ob- 
jects. RR Lyrae should be be identifiable on the basis 
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TABLE 1 

Properties of 838 LONEOS - I RR Lyrae Sample 



RA 


DEC 


< M L ON EOS V > 


Period 


Amplitude 


231.923019 


-22.413485 


14.48 


0.606328 


0.890 


234.829575 


-22.006363 


16.45 


0.644697 


0.734 


303.111725 


-23.159113 


15.18 


0.608629 


0.812 


302.516479 


-22.789234 


15.62 


0.670238 


0.705 


259.791656 


45.140347 


15.25 


0.468290 


1.083 



Note. — This table list a few properties of the 838 LONEOS - I RR 
Lyrae sample. The coordinates (RA,DEC) arc given in decimal degrees 
in the epoch of J2000.0. The mean magnitudes (< Mlon EOS v >) 
are tied to the Vsdss system. The periods are given in days. Finally, 
the amplitudes are measured from minimum to maximum brightness as 
defined by the RR Lyrae template that best fit the lightcurve. [The 
complete version of this table is in the electronic edition of the Journal. 
The printed edition contains only a sample] 




Fig. 9. — These color-color diagrams are constructed from SDSS 
photometry. The contours indicate the location of the stellar lo- 
cus in each of the 2D slices through the SDSS color-color space. 
The rectangular a reas are the RR Lyrae color boundaries are from 
llvezic e t al. (2005]). The red points are LONEOS - 1 RR Lyrae (78) 
with SDSS DR3 photometry. There are only two objects which do 
not pass the g — r cut, but are consistent with RR Lyrae stars in 
the three other colors. Thus, virtually all of these objects have 
colors consistent with RR Lyrae stars. 



of Sloan Di gital Sky Survey (SPSS) ph otometry alone 
(|Krisciunasll200H iKrisciunas et alJll998ft . 78 of the LO- 
NEOS - I R R Lyrae have clean S PSS Pata Release 3 
photometry (|Abazaiian et al.ll2005f ). Figure [5] shows the 
color-color diagrams constructed from SDSS PR3 pho- 
tometry. RR Lyrae stars are contained in a small region 
of SPSS color-color space (jlvezic et al.ll2005l ) as shown by 
the rectangular areas outlined in the figure. The contours 
indicate the location of the stellar locus in each of the 
2P slices through the SPSS color-color space. The red 
points are LONEOS - I RR Lyrae stars (78) with SPSS 
PR3 photometry. There are only two objects which do 
not pass the g — r cut, but are consistent with RR Lyrae 
stars in the three other colors. Thus, virtually all of these 
objects have colors consistent with RR Lyrae stars. We 
conclude that our RR Lyrae sample has a low degree of 
contamination. 

The survey contains fields with Galactic latitude as 
low as 10 degrees. Thus, there is a possibility that RR 
Lyrae from the thick disk have contaminated our sam- 
ple. Our period-amplitude cut likely has removed most 
of the thick disk RR Lyrae, but let us investigate how 
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Fig. 10. — This figure shows the predicted fraction of thick disk 
stars relative to halo stars which is expected in the LONEOS - I 
survey volume as a function of Galactocentric Radius. 

much thick disk contamination our sample might have 
using Galactic models. We use a double exponential to 
model the thick disk wi th a scale hei ght of 0.7 kpc and 
scale length of 3.0 kpc (|Lavdenlll995f ). These scales are 
based on measurements using RR Lyrae stars. Our halo 
is modeled as a power law with an index of —2.43, which 
is the value we found using the LONEOS - I RR Lyrae 
(See the next section for details). We assume that the 
thick disk is 10 times more populous than the stellar 
halo in the Solar Neighborho od, and the measured stel- 
lar number density is O.lpc -3 (|Reid fc Haw lcv 2000(1 . We 
also assume that the frequencies of RR Lyrae stars in the 
thick disk and halo are identical. This is a conservative 
assumption since thick disk RR Lyrae stars frequencies 
will be suppressed since thick disk RR Lyrae stars are 
more metal-rich. Given these parameters, we calculate 
the anticipated fraction of thick disk stars relative to halo 
stars in the LONEOS - 1 survey volume. Figure ITD1 shows 
this fraction as a function of Galactocentric Radius. This 
figure shows that the fraction of thick disk stars to halo 
stars is about 12% even at the smallest radii, where the 
largest contamination is expected. This low fraction of 
thick disk stars and our period-amplitude cut give us 
confidence that the contamination of thick disk stars in 
our sample is negligible. 

In order to investigate how well we recover various RR 
Lyrae light curve parameters, we simulated RR Lyrae 
light curves with known parameters and examined the re- 
covered parameters using the RR Lyrae period and light 
curve fitting procedure described above. We used the 
actual cadence for each subfield to generate simulated 
time series of RR Lyrae stars using template RR Lyrae 
light curves from Layden (1998) with various mean mag- 
nitudes. We used a know n RR Lyrae period- amplitude 
distribution (|Smithlll995f ) to generate the simulated time 
series with random phase. We assigned photometric er- 
ror bars and smear the simulated light curves based on 
the magnitude dependence of the LONEOS - I photo- 
metric errors. In Figure 111! the residuals of the input 
and recovered periods are shown. Periods are recovered 
with good accuracy. The fraction of simulated RR Lyrae 
contained within ±0.01 days are: 0.92, 0.90, 0.81, 0.70, 
for mean magnitudes of 15, 16, 17, and 17.5, respectively. 
We will investigate the possibility of misidentified peri- 
ods in the next section when we discuss the Oosterhoff 
components. 

Similarly, Figure [12] shows the residuals of the input 
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Fig. 11. — We have simulated RR Lyrae light curves with a range 
of periods and examined the periods recovered using our phasing 
and light curving fitting described in Section 3.1. The histograms 
show the difference between the input and recovered periods for 
four mean magnitudes. The fraction contained within ±0.01 days 
are: 0.92, 0.90, 0.81, 0.70, for mean magnitudes of 15,16,17, and 
17.5, respectively. 




a (Arrp) a (Amp) 

Fig. 12. — These histograms show the the residual between the 
input and recovered amplitudes from simulated RR Lyrae light 
curves of various mean magnitudes. The residual increases with 
mean magnitude. In addition, the distributions are skewed; the 
recovered amplitudes tend to be underestimated. 

and the recovered amplitude. The residual increases 
with mean magnitude. In addition, the distributions are 
skewed; the recovered amplitudes tend to be underesti- 
mated. This is a consequence of the asymmetric shape 
of ab-type RR Lyrae light curves, and the sparseness of 
our sampling. However, we find that the underestima- 
tion of amplitudes does not depend on the input period. 
Thus, underestimating the amplitude will only have the 
effect of shifting RR Lyrae stars downward in a period- 
amplitude diagram. We will discuss this more in the next 
section. 

3.2. Completeness 

In this section, we estimate our efficiency in detecting 
RR Lyrae stars (i.e., the completeness of the sample). 
The completeness depends on a number of RR Lyrae 
parameters, but most strongly on their magnitude and 
amplitude. One of our scientific goals is to measure the 
number density of RR Lyrae stars as a function of posi- 
tion in the Galaxy. This requires that we estimate the 
detection efficiency as a function of Galactic position. 

We divided our completeness analysis into photometric 
and temporal parts. In order to estimate the photomet- 
ric efficiency, we take a clean, unblended, and dereddened 



>, 1.0 




SDSS 



Fig. 13. — We have taken a clean sample of unblended SDSS DR3 
stars and queried the LONEOS - I database for a detection in each 
image in the database. This figure shows the overall photometric 
completeness (image-by-image) assuming SDSS is complete. Note 
that non-detection also include measurements that were flagged as 
being problematic. 

sample of stellar objects from SDSS DR3, which are con- 
tained in the LONEOS - I footprint. Since the SDSS 
detection limits are at least 2 magnitudes deeper than 
LONEOS, we can use SDSS to understand the degree to 
which LONEOS - I is photometrically complete. Using 
this sample, we performed a proximity query using our 
LONEOS - I database. For each SDSS star, we search 
in the LONEOS - I database for objects within 1.8 arc- 
sees. A total of 58,591 SDSS stars are queried in the 
LONEOS - I database. Next, we count the number of 
detections and non-detections of the object in every LO- 
NEOS - I image. A measurement is considered a non- 
detection if any one of the following criteria are true: 

• Signal-to-noise ratio is below threshold (5cr) 

• Measurements are flagged by photometric pipeline 
as bad (e.g., blended, pixel defects, non-stellar 
shape, cosmic ray hit) 

• More than one LONEOS - I object is detected, 
which indicates problems with time series gener- 
ation. 

These criteria are also used for the later analysis. 

Thus, this comparison with SDSS not only gives us 
an estimate of the photometric completeness based on 
signal-to-noise, but also an estimate of the efficiency of 
the photometric pipeline. The last two non-detection 
criteria are generally independent of magnitude. Thus, 
these will contribute to an overall reduction across all 
magnitudes. This is why the efficiency converges to only 
~ 90% in the bright magnitude range. Figure [T3l shows 
the LONEOS - I photometric completeness. We use 
the following parameterization for photometric efficiency 
shown in Figure [T51 

KVSDSS) = y„L-i8.™ - - U - ( 19 ) 

1 + e 0.47 

In addition, we did not find a significant dependence on 
the photometric completeness curves with position on 
the sky. Thus, we use this as our universal photometric 
completeness function. 

Using this parameterized photometric completeness, 
we performed the temporal completeness analysis. We 
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Fig. 14. — The RR Lyrae star completeness (photometric and 
temporal) for four fields with different number of images are shown. 
Note that even though fields 1354B and 1374B have identical num- 
bers of images, the completeness curves are different. This is caused 
by the different cadences of these two fields. These curves are used 
to correct the effective volume in the density measurements. 

generated simulated RR Lyrae times series in the same 
manner as described in Section 3.1, except that we in- 
corporated the parameterized photometric completeness, 
which provided a detection probability as function of 
magnitude. We randomly rejected individual time se- 
ries measurements based on this detection probability. 
Given the large number of subfields (~ 15, 000), we per- 
formed this analysis on one "typical" subfield per field, 
which had the median number of images; the number of 
images for subfields in a given field does not vary much. 
Thus, we perform this analysis on 421 "typical" subfields. 
We generated 1000 simulated time series for each "typi- 
cal" subfield and for each magnitude. All the simulated 
RR Lyrae times series were analyzed using the same se- 
lection criteria (see Figure SJ) as those used to generate 
the real LONEOS - I RR Lyrae catalog. The complete- 
ness functions were generated based on the percentage 
of simulated time series which passed our RR Lyrae se- 
lection criteria. The completeness functions as a func- 
tion of magnitude (M) for the "typical" subfields were 
parametrized by the following function: 

f'(M LO NEOS v) = TTZ^^TVIF ~ Z > ( 20 ) 

1 + e t 

where /i is the characteristic magnitude cutoff, T is the 
width of the cutoff, and Z is a horizontal offset. FigurefT4l 
shows the detection efficiency as a function of magnitude 
for four fields. All magni tudes have been corr ected for 
Galactic extinction using ISchleeel et all (|1998l ). These 
efficiencies are used when we calculate the survey vol- 
ume and corrected RR Lyrae star number densities in 
the following section. The completeness function for the 
"typical" subfield is used for all subfields of the same 
field. A histogram of the number of images in the "typ- 
ical" subfields is shown in the top left plot of Figure [T5l 
The median number of images is 31. In order to give a 
sense of the dependence of the completeness on the num- 
ber of images, Figure [15] shows the correlations of the 
magnitudes at 50% completeness and the completeness 
at the bright end (1-Z in Equation [20]) with the number 
of images in the top right and bottom left plots, respec- 
tively. The completeness is generally correlated with the 
number of images. 

4. ANALYSIS 



Fig. 15. — This figure shows various parameters and correla- 
tions related to the RR Lyrae star completeness (photometric and 
temporal). First, the top left plot is a histogram of the number of 
images for the typical subfields. The median number of images is 
31. The top right scatter plot shows the d epen dence of the mag- 
nitude at the bright end (1-Z in Equation I20| l on the number of 
images. This gives of sense of the completeness for bright stars, 
where signal-to-noise ratio is high. On the other hand, the bot- 
tom left scatter plot shows the dependence of the magnitude at 
50% completeness on the number of images, which gives a sense 
of the completeness towards the faint end. Finally, the bottom 
right histogram shows the distribution of the magnitude at 50% 
completeness. 

In this section, we explore the bulk properties of the 
LONEOS - 1 RR Lyrae star sample that was described in 
the previous section. First, we present the spatial distri- 
bution of the entire sample. Next, we present the period- 
amplitude distribution of the sample. Finally, we explore 
the spatial distributions of the different populations that 
emerge from the period-amplitude distribution. 

4.1. Radial Density Distribution 

As described in previous sections, we tie our abso- 
lute photometry to synthetic SDSS-based magnitudes 
(Vsdss)i Mloneos v- This allows for a straightforward 
calculation of distances from the intensity-weighted mean 
magnitudes. Distances are calculated from the following 
formula: 



We will use RR Lyrae absolute m agnitudes of My 
■ 0.71 ± 0.12 (jLavden et alJ I1996T ) derived from statis- 
tical parallax analysis for a sample with ([Fe/H]) = 
— 1.61. Assumin g that the Sun is at —7.8 kpc 
(|Carnev et al.l[l995] ) from the galactic center along the y- 
axis (xgauVgah z gai = (0,-7.8,0) kpc), we calculate the 
galactocentric 3-D coordinates of each RR Lyrae star. 

We measure the spherically symmetric radial density 
distribution by counting the number of RR Lyrae stars 
in galactocentric radial bins and calculating the survey 
volume contained in each bin. We use 1 kpc bins for 
all the density distributions. This provides us with a 
measure of the radial density distribution averaged along 
the survey lines-of-sight. Using the boundaries of the 
survey we numerically integrate the volume contained in 
each radial bin, taking into account detection efficiency 
at each position in the Monte Carlo integration. Figure 
[in] shows the survey effective volume distribution. 

Assuming a spherically symmetric distribution, Figure 
[l7]shows the radial profile for our sample with mean mag- 
nitudes in the range [13,17.5] with M v = +0.71. The 
one-sigma, asymmetric Poisson error bars are assigned 
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Fig. 16. — This figure shows the raw (circles) and the detec- 
tion efficiency-corrected (triangles) survey volumes. Note that the 
difference between the two curves is not monotonically increas- 
ing because there are lines-of-sight which are pointed towards the 
Galactic Center and also the anti-Center which contribute at dif- 
ferent ranges of Galactocentric radius. If the volume calculation 
assumes that the Sun is at the center of the Galaxy, the difference 
between these two curves does indeed increase monotonically. 




Fig. 17. — The radial density distribution of 838 afe-type RR 
Lyrae with a mean magnitude range from 13.0 to 17.5, using My = 
+0.71. We do not find any statistically significant deviations from 
a smooth stellar halo. Our survey volume does not intersect known 
locations of the tidal debris from the Sagittarius dwarf galaxy. 



and calculated using; the method of lGehrelsl (fl986h . Us- 
ing non-linear least squares, we fit power laws to the 
radial density distributions: 



p(R ga i) = (397 ± 56) x R { -J A3±om) (kpc 



(22) 



with x 2 per degree of freedom of 0.91 . All fitted param- 
eter estimates ar e calculated using a Monte Carlo boot- 
strap technique (|Press et al.lll992f i. Tabled] summaries 
results from other surveys. Thi s result is within 1 — a 
of the QUEST-I (jVivas fc Zinnl ( jQOl), M v = +0.55 , 
R SUN = 8.0 kpc) and the SDSS (jlvezic et all (|2000h . 
My = +0.7, Rsun = 8.0 kpc) results, assuming a 
spherical halo distribution. The power law exponent is 
slightly sensitive to the distance to the Galactic Cen- 
ter. If Rsun = 8.5 kpc, then exponent is reduced to 
—2.50 ± 0.05. As shown in Table [2 the measured dis- 
tance to the Galactic Center has been decreasing with 
time. This means that earlier measurements will tend to 
have steeper profiles than our measurement. 

The shapes of the stellar and dark matter halos should 
have encoded in them some information of their forma- 
tion process. In general, the isodensity contours for the 
RR Lyrae number density can be described by an ellip- 






Fig. 18. — Radial density profiles for spherical and variable flat- 
tening models are shown in the left an d right columns , respectively. 
The variable flattening model is that of Preston et al. ( 1991), where 
the flattening ratio variable from 0.5 to 1.0 with increasing distance 
from the Galactic Center. The top row of plots are for the entire 
sample of RR Lyrae stars. The second and third rows are for sub- 
samples of RR Lyrae, classified as Oosterhoff I and II, respectively. 
The LONEOS - I RR Lyrae data are fit by the spherical and vari- 
able flattening models equally well. The fits are not statistically 
different, 
soid: 



b 2 



1. 



(23) 



Let us consider the case where a 
and define the flattening ratio as 
isodensity contours can be described as: 

z 2 

x 2 + V 2 + ~o = a 2 , 



= b (i.e., a spheroid) 
q = c/a. Thus, the 



(24) 



where q < 1 and q > 1 are oblate and prolate spheroids, 
respectively. It is widely observed that the inner regions 
{Rgai — 5 — lOkpc) of the stellar halo are highly flat- 
tened (q ~ 0.5), while the outer regions are more spheri- 
cally symmetric (e.g., [Kinman et al.|[l966T : iPreston et al.l 
[19911: iChiba k Beers! 12000 ) . We fit spheroids with var- 
ious constant flattening to the LONEOS - I RR Lyrae. 
In addition, we also use a model wit h a spatial variable 
flattening ratio (|Preston et al.|[l99lT ). where q varies as 
function of semi-major axis, a: 



q = 



_ / 9o + [1 - qo]jo if a < 20 kpc 
1 if a > 20 kpc 



where qo = 0.5 and a is the semi-major axis. In this 
variable q model, the flattening smoothly varies from 0.5 
to 1.0 over the course of 20 kpc. Figure [T8l shows the 
radial profiles and power law fits for both the spherical 
and the variable flattening models. The power law ex- 
ponents, x 2 i degrees of freedom, and % 2 probabilities 
are listed in Table [3] The degrees of freedom are not 
constant since the data are binned and the value of a 
depends on q. All the flattening models have larger 
X 2 probability than the spherically symmetric model in- 
dicating that flattened models are favored, but there is 
no clear minimum in the x 2 probabilities. The variable 
q model yields a steeper density profile with a power law 
exponent of —3. 15 ±0.07 with a x 2 P er degree of freedom 
of 0.7. This agrees with most other RR Lyrae surveys 
(Table 

The poor constraint on the stellar halo flattening from 
the LONEOS - I RR Lyrae is likely a result of the fact 
that the Galactic latitude distribution of these RR Lyrae 
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Halo Hit Lyrae Survey Comparisons 




Survey 


Tracer 


Exponent 


Flattening Ratio 


<q) 


Objects 


Ro (kpc) 


Comment 


LONEOS - 


I (All) 


RR Lyrae 


-2.43 ± 0.06 


q = 1 




838 


7.6 


M v = ±0.71 


LONEOS - 


I (All) 


RR Lyrae 


-3.15 ± 0.07 


Variable q 




838 


7.6 


M v = ±0-71 


LONEOS - 


I (Oo I) 


RR Lyrae 


-2.26 ± 0.07 


q = 1 




619 


7.6 


M v = ±0.71 


LONEOS - 


I (Oo II) 


RR Lyrae 


-2.88 ± 0.11 


q = 1 




219 


7.6 


My = ±0.71 


QUEST-I ( 


Vivas & Zinn (2006)) 


RR Lyrae 


-2.5 ±0.1 


q = 1 




395 


8.0 


My = ±0.55 


QUEST-I ( 




RR Lyrae 


-3.1 ± 0.1 


Variable q 




395 


8.0 


M v = ±0.55 


SDSS flvez 




RR Lyrae 


-2.7 ± 0.2 


q = 1 




148 


8.0 


M v = ±0.7 


Chiba & tJ 


eers (2000) 


Halo «[Fe/ff]> < -1.8) 


-3.55 ± 0.13 


Variable q (qq ~ 


0.65) 


413 




Kinematics 




sers (2000) 


Halo (-1.6 < ([fe/H]) < -1.0) 


-3.47 ± 0.18 


Variable q (qq ~ 


0.55) 


331 




Kinematics 


Wetterer & 




RR Lyrae 


-3.0 ± 0.08 


q = 1 




42 


7.6 


M v = ±0-74 


Wetterer & 


McUraw (TSSo - ) 


RR Lyrae 


-3.5 ± 0.08 


Variable q 




42 


7.6 


M v = ±0.74 


Preston et 


al. (1991) 


RR Lyrae 


-3.2 ±0.1 


Variable q 




47 


8.0 


M v = ±0.6 


Saha (1985 




RR Lyrae 
Globular Clusters 


-3.0 
-3.5 


q = 1 




29 
121 


8.7 


P-L-A relation (Mg = M bot ± 0.34) 


^a^m^a984) 


RR Lyrae 


-3.1 ± 0.2 


q < 0.9 




24 


8.7 


M B = ±0.9 


Oort & Plaut (1975) 


RR Lyrae 


-3.0 


q < 0.8 




1108 


8.7 


M B = ±0.9 


Kinman ct 


al. (1966) 


RR Lyrae 


-3.5 


q= 0.6 - 0.S 


1 


38 


10.0 


M v = ±0.42 



Note. — We < 



the LONEOS - I RR Lyra 



TABLE 3 
Flattening Fits 



Sample 


Flattening 


Exponent 


x 2 


Degrees of 


Probability 




Ratio (q) 






Freedom (v) 


(Q(x 2 ,■>)) 


All 


1.0 


-2.43 


24.7 


24 


0.42 


All 


Variable q 


-3.15 


18.3 


23 


0.74 


All 


0.9 


-2.48 


16.1 


24 


0.88 


All 


0.8 


-2.53 


17.7 


24 


0.82 


All 


0.7 


-2.57 


14.4 


23 


0.91 


All 


0.6 


-2.58 


19.8 


23 


0.65 


All 


0.5 


-2.57 


13.4 


22 


0.92 


All 


0.4 


-2.57 


15.4 


22 


0.84 


All 


1.1 


-2.35 


17.5 


25 


0.86 


All 


1.2 


-2.29 


19.1 


25 


0.79 


Ool 


1.0 


-2.26 


18.9 


24 


0.76 


Ool 


Variable q 


-2.93 


16.0 


23 


0.86 


Ool 


0.9 


-2.29 


13.8 


24 


0.95 


Ool 


0.8 


-2.36 


14.8 


24 


0.93 


Ool 


0.7 


-2.40 


12.6 


23 


0.90 


Ool 


0.6 


-2.42 


19.8 


23 


0.65 


Ool 


0.5 


-2.41 


10.5 


22 


0.98 


Ool 


0.4 


-2.41 


10.9 


22 


0.98 


Ool 


1.1 


-2.18 


13.3 


25 


0.97 


Ool 


1.2 


-2.15 


15.9 


25 


0.92 


OoII 


1.0 


-2.88 


7.2 


24 


0.9996 


OoII 


Variable q 


-3.70 


5.9 


23 


0.9998 


OoII 


0.9 


-2.97 


5.6 


24 


0.9999 


OoII 


0.8 


-3.02 


5.7 


24 


0.9999 


OoII 


0.7 


-3.09 


7.5 


23 


0.9990 


OoII 


0.6 


-3.07 


5.8 


23 


0.9998 


OoII 


0.5 


-3.01 


9.6 


22 


0.9895 


OoII 


0.4 


-3.04 


13.0 


22 


0.9331 


OoII 


1.1 


-2.80 


10.0 


25 


0.9966 


OoII 


1.2 


-2.71 


11.8 


25 


0.9881 


Note. — 


Power law 


fits with vary: 


mg flattc 


ning models. 





are highly concentrated at mid-range Galactic latitude 
with no fields near the Galactic Poles. The difference 
between the spherical and flattened model is smallest at 
lower Galactic latitudes. 

All the recent RR Ly r ae su rveys are in disagree- 
ment with IChiba fc Beersl ((2000) who used ~ 1200 non- 
kinematically selected metal poor stars. One reason for 
the difference could be the method used. RR Lyrae sur- 
veys use 3D Galactocentric coordinates to measure den- 
sities. In these surveys, the completeness and survey vol- 
ume a re relatively well-defined. However. lChiba fc Beersl 
(2000) use data derived from many different surveys and 
is not a volume-limited sample. They use a maximum 
likelihood method which relies on kinematic information 
and Jeans' theorem to measure the power law. Since 
the completeness and survey volumes are not easily un- 
derstood, it is not possible to measure the density dis- 
tribution using a binning technique. Thus, it is diffi- 
cult to understand if there is a real difference between 
RR Lyrae stars and generic metal-poor stars. In addi- 
tion, RR Lyrae stars and the metal poor halo stars from 
I Chiba fc Beersl (|2000l) may possibly be tracing different 



surveys using RR Lyrae stars and other tracers. 

populations. RR Lyrae abundances are influenced by the 
horizontal branch morphology of the RR Lyrae progeni- 
tor system. Thus, they may not reveal the entire history 
of Galactic formation and could have different spatial 
distributions than a sample of metal poor halo stars. 

4.2. Period- Amplitude Distribution 

The period-amplitude distribution of RR Lyrae stars 
can be used to explore the for mation history of the Galac- 
tic halo (e.g., |Catelan 2006). If some fraction or all of 
the stellar halo was formed from the accretion of dwarf 
galaxies, which might have very diverse chemical histo- 
ries and ages (i.e., different horizontal branch morpholo- 
gies), the signatures of these accretion events should be 
imprinted in the period-amplitude distribution. Thus, 
accretion activity can be probed even after phase-space 
has been sufficiently mixed that spatial coherence is lost. 

Our large sample of RR Lyrae stars allows us to fur- 
ther investigate this possibility and to explore differences 
in the two components. In lower-right panel of Figure 
[7J we plot the period-amplitude distribution of 838 ab- 
type RR Lyrae stars. Since LONEOS uses an unfiltered 
CCD system, whose sensitivity peaks in the red, ampli- 
tudes for our RR Lyrae will be systematically smaller 
than those taken in the V-band. The observed ampli- 
tude of RR Lyrae increases as you observe them in bluer 
passbands. 1 Using a subset of LONEOS - I RR Lyrae 
(12 stars) which have previously measured V-band am- 
plitudes taken from SIMBAD database, we have empiri- 
cally derived a scaling relation between the LONEOS am- 
plitudes (Aloneos) and V-band amplitudes (A v ) and 
scaled the M3 RR Lyrae period-amplitude relation to 
the LONEOS - I amplitude scale. The mean value of 
Aloneos/A v is 0.7 with a standard deviation of 0.1. 
Figure [19] shows that data from the lower-right panel in 
Figure [7j overlaid with the linear period-amplitude re- 
lation for M3 (the q uintessential Ool globular cluster) 
(|Cacciari et al.ll2005[ ) corrected to the LONEOS - I am- 
plitude scale. Despite only having 12 stars to scale LO- 
NEOS amplitudes to V-band amplitudes, the main locus 
agrees well with the period-amplitude relation of M3. 
Thus, we associate the main locus with the Ool compo- 
nent of the stellar halo. In addition, Figure [TU shows 
a striking concentration of RR Lyrae (26% of the total 

1 Note that the transformation to the V$DSS magnitude sys- 
tem does not affect the measured LONEOS - I amplitudes. The 
amplitudes are strictly a function of the passband used. 
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Fig. 19. — The LONEOS RR Lyrae period-amplitude distribu- 
tion (same as Figure[7} with the linear log(period)-amplitude rela- 
tion for globular cluster M3 (Ool) superimposed (solid line). Note 
the amplitudes are based on the unfiltered response of the LO- 
NEOS - I camera system. The M3 relation has been converted to 
the LONEOS - 1 amplitude scale. We derived a conv ersion between 
LONE OS - 1 amplitude and V-band amplitudes from lCacciari et al.l 
(2005). The dashed line indicates where we divide our sample into 
Ool and OoII components. 




Fig. 20. — Histogram of the period shift of RR Lyrae relative 
to the primary locus (Ool). We fit a line to the primary locus and 
measure A Log(Period) which refers to the Log(Period) distance 
(at fixed amplitude) of each RR Lyrae star from this best fit line. 
Thus, the primary locus is centered at A Log(Period) = 0.0. The 
location of the secondary locus is located at a pproximately the lo- 
cation where we wou ld expect the OoII locus (ISandage et al.|[1981l ; 
ICacciar i et al. 2005]) . For the subsequent analysis, we divide our 
sample at A Log(Period) = 0.045, where period is measured in 
days. 




-0.10 -0.05 0.00 0.05 0.10 0.15 0.20 

A Log(Period) 

Fig. 21.— Period shifts (with respect to M3) of the QUEST-I RR 
Lyrae Survey is shown as a solid line, while the LONEOS-I period 
shifts are shown as a dashed line. In contrast to the LONEOS- 
I sample, there is no clearly displaced secondary peak associated 
with an OoII component for the QUEST-I sample. However, there 
is a long period tail. 

ponent, and the shifted locus with an OoII halo compo- 
nent. 

In Figure [21] we plot the period shift dis tribution of 
395 R R Lyrae from the QUEST-I catalog (| Vivas et all 
200 1 ■ with mean magnitudes in the range [13.5, 19.7] as 
well as the LONEOS-I period-shift distribution. This 
QUEST-I distribution also contains a sizable fraction 
(~ 20 %) of longer pe riod RR Lyrae than the main Ool 
locus. Catelan (2006|) has pointed out that the Ooster- 
hoff classification the QUEST-I RR Lyrae may be unre- 
liable due to their sparse temporal sampling. This could 
cause a smearing of the period shift distribution in the 
QUEST-I data. 

It has been well established that Ool and OoII globu- 
lar clusters show spatial and kinematic differences, which 
may signify different origins (|Lee fc Carnevl [l999). It is 
natural to ask if the Ool and OoII components of the stel- 
lar halo show similar differences. We divided our sample 
into the two components described above. The division 
was positioned at the minimum of the period shift dis- 
tribution, A(P) = 0.045. The division is denoted by the 
dashed line in Figure [THJ The radial number density pro- 
files for each component are shown in Figure [H] From 
these profiles, we derive the following radial dependences: 

Pool ~ iT 2 - 26±0 - 07 , (26) 



sample) with longer periods offset from the main locus 
of points. 

We fit a line to the primary locus (P pr imary{A)) and 
measure the period shift, 



A Log(Period) = P RR (A) - P 1 



primary 



(A), (25) 

which refers to the Log(Period) distance (at fixed ampli- 
tude) of each RR Lyra (Prr(A)) from this best fit line 
(A = —4.76 x LogiPprimary) — 0.51). Figure [20l shows the 
distribution of the period shift at fixed amplitude from 
the main locus. This figure clearly shows the existence 
of a secondary locus. It is shifted from the main locus 
by ALog(P) ~ 0.075 2 , which is approximately the shift 
between RR Lyrae in Ool to OoII globular clu sters (e.g., 
ISandage et aTlll 9 811 : [clement fe Sheltonlll999f) . Thus, we 
associate the main locus of points with an Ool halo com- 

2 Periods are measured in days. 



n-2.88i0.ll / 07 \ 
POoII ~ ti . (Z{) 

The OoII component has a steepe r profile than the Oo l 
component. This is consistent with lLee fc Carnevl (|1999D 
who, using a smaller RR Lyrae sample, found that the 
Ool component dominated farther from the Galactic 
plane, possibly indicating that the OoII component falls 
off more rapidly than the Ool component. In the fol- 
lowing sections, we will examine whether this difference 
could be caused by some systematic effect in the detec- 
tion of the two populations. 

The possibility of an intermediate Oosterhoff halo com- 
ponent still exists, but our current sample cannot address 
this. It is critical to effectively identify and veto RR 
Lyrae which exhibit the Blazhko effect (i.e. amplitude 
and phase variations). These stars will add scatt er to 
the period-amplitude relation (Cac ciari et all 120051) . In 
addition, we have found that the LONEOS - I measured 
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Fig. 22. — The spherically symmetric (q=l) radial distribution 
for the Ool and OoII components. The bins with zero density are 
included in the f it; the l a Poi sson upper limits are calculated using 
the method of lGehrelsl il986ft . 

RR Lyrae amplitudes tend to be underestimated (Figure 
fl~2|) . This effect would tend to diminish the number of 
RR Lyrae with intermediate Oo classification. 

Additionally, microlensing surveys (e.g., MACHO, 
OGLE) have found that as much as a third of RR Lyrae 
stars in the LMC and S MC exhibit the Blazhko effect 
(e.g.. lAlcock et al1 [2000). In order to address the possi- 
bilities of an intermedia te Oo compon ent, we need better 
temporal coverage (e.g. JCate lan 2006). A mixture of Ool 
and OoII could masquerade as an intermediate. 

Is the difference in radial distributions associated with 
the Oosterhoff effect? Does the Oosterhoff effect mani- 
fest itself in such a way? Previously, we have chosen to 
divide our full sample into two components that resem- 
ble Ool and OoII populations; this maximizes the role 
of the Oosterhoff effect on the sample. Instead, we now 
divide our entire sample along a line in period-amplitude 
space perpendicular to Oosterhoff dividing line, mini- 
mizing the role of the Oosterhoff effect. The resulting 
power law exponents (a) are consistent with each other 
(a = -2.50 ± 0.07 and a = -2.39 ± 0.1). This indi- 
cates that the Oosterhoff effect is likely associated with 
differing radial distributions. 

5. DISCUSSION 

As described in the previous section, the LONEOS - I 
RR Lyrae data have two interesting characteristics. The 
Oosterhoff effect is manifest in the Galactic stellar halo. 
Moreover, the two Oosterhoff components we have iden- 
tified have distinct radial distributions. 

5.1. Potential Systematic Errors and Biases 

In order to test whether the the Oosterhoff effect in 
the stellar halo and difference in the radial distribution 
of the two Oosterhoff components is real, we performed 
a series of tests to examine systematic effects that might 
account for the difference. It appears that this result is 
robust. 

Period Errors — Let us assume that the stellar halo con- 
sists of solely an Ool population of RR Lyrae stars. Is 
it possible that Ool RR Lyrae stars might be misiden- 
tified as OoII? For example, it is possible that periods 
are grossly mis-estimated. Could our timeseries anal- 
ysis and detection process artificially produce an OoII 
population? We simulated RR Lyrae light curves with 
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Fig. 23. — We simulated RR Lyrae light curves with period and 
amplitude distributions of Ool and OoII RR Lyrae stars. These 
simulated RR Lyrae light curves were processed through our de- 
tection pipeline. This figure shows the period shift histograms for 
one field with 28 images (1374B). We have scaled the OoII input 
sample to contain 26% of the Ool component, which is what we 
observe in the LONEOS - I sample. The detection process tends 
to broaden the period shift distributions and produces long tails. 
However, the Ool and OoII are still recognizable after the detection 
process. 

period and amplitude distributions of Ool and OoII RR 
Lyrae stars. These simulated RR Lyrae light curves were 
processed through our detection pipeline, which was de- 
scribed in the previous section, and used in the selection 
of final LONEOS - I RR Lyrae sample discussed in this 
section. Figure [23] shows the period shift histograms for 
one field with 28 images (1374B). We have scaled the 
OoII input sample to contain 26% of the Ool compo- 
nent, which is what we observe in the LONEOS - I sam- 
ple. The detection process tends to broaden the period 
shift distributions and produces long tails. However, the 
Ool and OoII components are still recognizable. 

Efficiency — In order to calculate the RR Lyrae star 
number density, we need to estimate the detection ef- 
ficiencies. This was described in the previous section. In 
calculating the detection efficiencies, we averaged over 
the period and amplitude distributions from previous 
surveys. These distributions are essentially for the Oost- 
erhoff I class. The assumption that we have made is 
that the detection efficiencies for Ool and OoII are iden- 
tical. In order to investigate differences in the detec- 
tion efficiencies, we simulate RR Lyrae light curves to 
explore the detection efficiencies in the period-amplitude 
space. We performed more extensive simulations exam- 
ining only fields with less than 30 images, which are fields 
with the fewest images and represent the worst-case sce- 
nario. We simulated 310,000 RR Lyrae stars with mean 
magnitudes of 15, 16, 17, & 17.5 (Figure EJ) . There is 
a strong trend with amplitude which is expected given 
the photometric error trend with magnitude. We also 
observe a deficiency near a period of 0.5 days (-0.3 on 
the log scale), which is caused by the sampling window 
function of these fields. We find no dependence of the 
efficieny with period at fixed amplitude. Thus, we are 
confident that the efficiency differences do not account 
for the difference in the the radial distributions of the 
Ool and OoII components. 

What effect could an error in the detection efficiency 
have? We have measured the power law indices for two 
Oosterhoff components without including the detection 
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Fig. 24. — The detection efficiency in period-amplitude space 
for fields with less than 30 images, which are fields with the fewest 
images. We simulated 310,000 RR Lyrae stars with mean magni- 
tude of 17.0. Note that color table has been inverted (i.e., black 
corresponds to 100% efficient). We did not find any dependence of 
the efficieny with period at a fixed amplitude. Also, one can see a 
deficiency at a period of 0.5 days (-0.3 on the log scale), which is 
caused by the window function of these fields. 





Fig. 25. — 2MASS colors for the Ool and OoII components. 
There is no significant color difference between the two compo- 
nents. 



efficiency correction. This yields the following fits: 



POol 



POoII 



R 



-2.68±0.06 



R 



-3.17±0.1 



(28) 



(29) 



The detection efficiency correction tends to make the dis- 
tributions shallower. Even if we take the extreme sce- 
nario where our detection efficiency correction is only ap- 
plied to the OoII component and we assume that the Ool 
component has a 100% detection efficiency, the difference 
in the radial distributions of the Oosterhoff components 
remains statistically significant. 

Color Differences — Is it possible to explain the Ool and 
OoII difference because of color difference? Does the 
photometric calibration introduce a bias based on color? 
We looked to see if our Ool and OoII components have 
different color distributions. First, we have acquired 
dereddened 2MASS colors for several hundred of the LO- 
NEOS - I RR Lyrae stars. Figure [25] show the color 
distributions for each Oosterhoff component. The dis- 
tributions appear to be consistent and no infrared color 
difference is seen. In addition, Figure [26] shows the color 
distributions for the 78 RR Lyrae stars with SDSS DR3 
photometry. While this sample is small, we do not see 
major differences in the SDSS colors of the LONEOS - I 
RR Lyrae stars, except for the u — g color, where the 
Ool component appears to be bluer. This u — g color 
difference could be due to a difference in metallicity. 

Zero-point Differences — Zeropoint differences between 
the two Oosterhoff components may arise in several 
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Fig. 26.— SDSS color distributions for LONEOS - I RR Lyrae 
stars (78). The Oo I component is shown as the solid curve, while 
the Oo II component is shown as a dashed line. 
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Fig. 27. — The dependence of the power law exponent on the 
RR Lyrae absolute magnitude for the entire LONEOS-I RR Lyrae 
sample, Ool and OoII components is shown. If we fix the absolute 
magnitude of the Ool stars to My = 0.71, this would require that 
the absolute magnitude of the OoII component be roughly My > 
1.5. This value is not within the observed range of My for RR 
Lyrae of any kind. Thus, an absolute magnitude difference between 
the Ool and OoII components cannot explain the different radial 
distributions. 

forms: absolute magnitude differences, photometric ze- 
ropoint offsets, and mean magnitudes errors. First, we 
explore whether the difference in the radial distribution 
of the Ool and the OoII components could result from 
intrinsically different absolute magnitudes. In Figure [2~T1 
we show the dependence of the power law exponent on 
the absolute magnitude for the two Oosterhoff compo- 
nents. If we fix the absolute magnitude of the Ool stars 
to My = 0.71, this would require that the absolute mag- 
nitude of the OoII component be roughly My > 1.5 in 
order to have the same power law slope. This value is 
not within the observed range of My for RR Lyrae of 
any kind. Thus, an absolute magnitude difference be- 
tween the Ool and OoII components cannot explain the 
different radial distributions. 

Does the photometric calibration introduce a bias 
based on color? Even though we do not see a clear color 
difference between the Oosterhoff components, let us ex- 
plore how a difference could affect our results. As we 
saw in Section 2 (Figure [2J) , the offset from LONEOS - I 
magnitudes (tied to GSC2 F) to Vsdss is a function of 
color. A systematic error in this offset directly trans- 
lates to a systematic offset in the mean magnitudes (and 
distances) of all the RR Lyrae stars. This will lead to 
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different radial profiles. We have used an offset for a 
color range of RR Lyrae stars. Figure [2] indicates that 
the offset may vary by ~ 0.3 at the extremes of the RR 
Lyrae color range. To give us an idea of how this affects 
the radial profiles of the Oosterhoff components, let us 
look at Figure l27l where we explore how changes in the 
absolute magnitude affect the radial distributions. An 
overall shift in the LONEOS - I zeropoint is essentially 
changing the RR Lyrae absolute magnitude. If we fix 
My and assume the two Oosterhoff components have a 
0.2 mag overall zeropoint difference, this can not explain 
the difference in the two radial distributions. One needs 
an overall zeropoint difference of ~ 1 mags, which we are 
confident does not exist. 

It is possible that intrinsic light curve shape differences 
in the Ool and OoII components could result in a sys- 
tematic bias in the mean magnitudes of the Ool and OoII 
components. Any systematic difference in the measured 
mean magnitude is a zeropoint offset. However, Figure 
1271 shows that the mean magnitude bias would have to 
be a sizable fraction of the RR Lyrae amplitude and thus 
unrealistic that the radial profiles could result from a bias 
in the measure of the mean magnitudes. 




Fig. 28.— We have divided the LONEOS - I RR Lyrae sample 
into regions above and below the Galactic plane. This figure shows 
the spherically symmetric radial distributions for the entire sample 
as well as for the Ool and OoII components. The OoII component's 
profiles are consistent with the full OoII sample. However, the 
Ool shows a 2.8 a difference above and below the Galactic plane 
with power law exponents of —2.61 ± 0.09 and —2.16 ± 0.13. If 
this difference is real, it could indicate spatial asymmetries in the 
Ool component. However, the above and below plane difference in 
power law slope is of marginal statistical significance. 



Spatial Sub-Samples — We have divided the LONEOS - 
I RR Lyrae sample into regions above and below the 
Galactic plane. This division roughly divides the entire 
sample in half (372 and 466, above and below the plane, 
respectively) . Figure [28] shows the spherically symmet- 
ric radial distributions for the entire sample as well as 
for the Ool and OoII components (see also Table |4|). 
The OoII component's profiles are consistent with the 
full OoII sample. However, the Ool shows a 2.8 a differ- 
ence above and below the Galactic plane with power law 
exponents of —2.61 ±0.09 and —2.16 ±0.13, respectively. 
If this difference is real, it could indicate spatial asym- 
metries in the Ool component. However, the above and 
below plane difference in power law slope is of marginal 
statistical significance. We are engaged in the analysis 
of the LONEOS II data set which (with its much greater 
sky coverage) should lay this issue to rest one way or 
another. It is also possible that this difference is caused 
by some systematic variation in the LONEOS - I pho- 
tometric completeness function which we have not taken 
into account. Recall, that we use a global photometric 
completeness function for LONEOS - I and a field-by- 
field temporal completeness correction. We have inves- 
tigated if there was a systematic difference between the 
LONEOS - I fields above and below the Galactic plane. 
First, we derived a photometric completeness function for 
the LONEOS - I fields above and below the plane using 
SDSS Data Release 5. Using the same criteria that we 
used to obtain the global photometric completeness func- 
tion (Figure [T3"]) . we queried the LONEOS - 1 database for 
a detection in each image in the database and computed 
photometric completeness functions above and below the 
Galactic Plane. There were about one million detections 
of SDSS DR5 stars in the LONEOS - I database above 
and below the Galactic Plane. However, we did not find a 
significant difference between photometric completeness 
functions above and below the Galactic Plane. Next, 
we performed lightcurve simulations, like those in Figure 
|2"41 to compare the temporal detection efficiency in the 
period-amplitude space for LONEOS - I fields above and 
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Power law fits above and below tbe Galactic Disk 







Sample (Size) 


Exponent 


*ga 


> 





All (372) 


-2.67 ± 0.11 


'go. 


< 





All (466) 


-2.43 ± 0.08 


All 






All (838) 


-2.43 ± 0.06 


"ga 


> 





Ool (270) 


-2.61 ± 0.09 


*ga 


< 





Ool (349) 


-2.16 ± 0.13 


All 






Ool (619) 


-2.26 ± 0.07 


*S» 


> 





OoII (102) 


-3.12 ± 0.19 


* ga 


< 





OoII (117) 


-3.09 ± 0-19 


All 






OoII (219) 


-2.88 ± 0-11 



below the Galactic plane. Again, no significant differ- 
ence was seen for the fields above and below the Galactic 
plane. However, the Ool and OoII still have different pro- 
files in each hemisphere, but the extent of the difference 
is not presently clear. 

CCD Focal Plane Distribution — In order to obtain the 
wide field of view needed for LONEOS, the telescope 
used has a Schmidt design. A Schmidt telescope con- 
sists of a spherical primary mirror and a corrector lens 
at its radius of curvature. However, distortions in the 
focal plane persist. These distortions increase near the 
edges of the focal plane (i.e., with increased distance from 
the optical axis). Such distortion might affect the am- 
plitude of our RR Lyrae stars. We examined the distri- 
bution of our RR Lyrae stars in the CCD focal plane to 
search for any unusual clustering that might result from 
the optical distortions. We use an extension of a one- 
dimensional Kolm ogorov- Smirnov test to two-dimensions 
(|Press et al.lll992D . We find that the RR Lyrae distri- 
bution is not significantly different from a random 2D 
distribution. In addition, we examine the distributions 
of the two Oosterhoff components. The KS probability 
indicates that these distributions are not significantly dif- 
ferent from one another. 

Large Amplitude Sample — The detection efficiency gener- 
ally decreases with decreasing amplitude. We examined 
the radial distributions of the Ool and OoII components 
while restricting the sample to include only RR Lyrae 
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stars with amplitudes greater than 0.6. The difference in 
the radial distributions of the two components remains 
unchanged. 

Increased Cadence Subsample — Finally, we look at a sub- 
sample selected from fields with more than 32 epochs. 
Such a subsample should be less affected by efficiency dif- 
ferences. This results in 370 Ool and 129 OoII RR Lyrae 
stars. The power law exponents are a = —2.46 ± 0.1 
and a = —2.95 ± 0.14 for the Ool and OoII components, 
respectively. Again, the difference in radial distributions 
remains. 

5.2. Implications for Galactic Formation 

The fir st serious model of galaxy formation was pro- 
posed by lEggen et al.l (|1962f ) (ELS). Based on correla- 
tions between metallicity, orbital eccentricity, angular 
momentum, and velocity perpendicular to the Galactic 
Plane in a sample of nearby, high velocity (halo) field 
stars, they concluded that the Galaxy formed from the 
rapid (within the free fall time) , monolithic collapse of an 
overdensity in the primordial density field. The first stars 
that formed during the collapse inherited the low metal- 
licity and predominantly radial motion of the collapsing 
gas, while conservation of angular momentum and energy 
dissipation led to the formation of a flattened and circu- 
larly rotating thin disk of gas. Metal-rich stars would 
form in this disk as the interstellar medium became en- 
riched with the products of nucleosynthesis during stellar 
evolution. 

There has been mounting evidence to suggest that 
the ELS scenario d oes not agree with observations. 
ISearle fc Zinnl (|1978l ) (SZ) presented an extensive survey 
of globular clusters in the outer halo that showed a spread 
in their horizontal branch (HB) morphology (i.e., color 
distribution of stars along the HB) that is independent of 
metallicity, the "second parameter problem" . This effect 
was not evident in clusters with smaller galactocentric 
distances. SZ proposed that age was the second parame- 
ter and that the outer clusters originated in satellite sys- 
tems that were subsequently accreted by the Milky Way. 
The spread in HB morphology can be explained by dif- 
ferent star formation histories in the accreted systems. 
Today, the most successful large-scale structure forma- 
tion models are based on hierarchical clustering (i.e., the 
Cold Dark Matter (CDM) paradigm). The SZ scenario 
for galaxy formation is attractive since it would fit within 
the CDM paradigm of structure formation. Several ac- 
creted system have been detected in the M ilky Way (e.g., 
Ilbata et all 120011 : iBelokurov et al J 120071 ). which would 
provide direct support for SZ scenario. 

In studies of globular cluster RR Lyrae stars, the lo- 
cation of an RR Lyra star in period-a mplitude space 
has been correlated t o metallicity (e.g.. ISandaed [19931 
IContreras et al.1 l2005f ) and age relati ve to the Zero- 
Age Horizontal Branch (Z AHB) (e.g.. Eee et al l [l990t 
Clement fc Sheltonl [l999h iClement fc Sheltonl fl999). 



Lee fc Carney! (|1999l ). and Uurcsik et al.l (|2003h conclude 



that the Oosterhoff effect is caused by evolution off 
the ZAHB. They found that OoII globular clusters arc 
not correlated with metallicity, but rather their RR 
Lyrae stars appear to have evolved off the ZAHB. Rela- 
tive age estimates between the O ol and OoII globular 
cluster populations are 2-3 Gyr (|Lee fc Carne 
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iLee fc Carnevl (|1999l) have found that the kinematic 
properties of Ool and OoII globular clusters are dif- 
ferent. They found that the Ool clusters have little 
or no net angular momentum, which points to a for- 
mation scenario like that of SZ. On the other hand, 
OoII clusters were found to have a net prograde rota- 
tion (Vrotation = —94 ± 47 km/s), which indicates that 
their formation scenario resembles the ELS scenario. The 
kinematic data and the relative age differences (indicated 
by the Oosterhoff type) of Ool and OoII clusters point to 
a dual mode formation scenario for the globular clusters. 

Our detection of the Oosterhoff effect in the Galactic 
stellar halo indicates that the formation of the globular 
clusters may be connected to that of the Galactic stellar 
halo. If the Oosterhoff type indicates an age difference, 
then the Oosterhoff effect in the stellar halo suggests a 
time gap between the formation of the Ool and OoII 
components. It should be noted that an age difference 
between Ool and OoII components does not give us any 
information about when the Ool component may have 
been accreted into the Milky Way. However, the pres- 
ence of the two distinct components in the stellar halo 
discovered in LONEOS - I indicates that the progeni- 
tors of the stellar halo had diverse ages and/or chemical 
properties. 

One way to connect the Ool and OoII components to 
the ELS and SZ scenarios would be to search for phase- 
space substructure. First, if the Ool component was 
formed predominantly in the SZ scenario than we would 
expect the stellar halo to have over- and under-densities. 
While there has been much excitement concerning the 
discovery of the overdensities in the stellar halo that are 
associated with accreted systems, the stellar halo should 
also have voids or underdensities if it was formed in an 
SZ scenario. The radial distributions of both components 
exhibit no significant substructure. However spatial sub- 
structure tends to be quickly erased. Thus, the lack of 
substructure in the Ool component may simply indicate 
that the number of accretions in the past 1-2 Gyrs has 
been low. In addition, LONEOS - I probes a small frac- 
tion of the Galactic Halo. A much larger survey is needed 
to explore whether there are differences in the amount of 
substructure in the Ool and OoII components. 

Besides the possible time gap between the formation of 
the Ool and OoII components, the difference in the ra- 
dial distributions of the Ool and the OoII components is 
suggestive of multiple physical mechanisms for the forma- 
tion of the stellar halo. If one assumes that the Galactic 
stellar halo was formed by a combination of ELS and SZ 
scenarios, then one might associate the two Oosterhoff 
components of the stellar halo to each of those scenar- 
ios. The relatively steeper radial distribution of the OoII 
component as compared to the Ool component could be 
a manifestation of the fact that an ELS scenario is more 
dissipative than an SZ scenario. An ELS scenario would 
be more dissipative since there is much gas in the system. 
Loss of energy (e.g., from radiative transfer phenomena) 
would result in a contraction of the halo. Whether the 
contraction was isotropic or not depends on the amount 
of angular momentum in the system. 

However, our detection of the Oosterhoff effect in the 
stellar halo does not necessarily imply that the stellar 
halo was formed through a combination of ELS and SZ 
scenarios. It might be able to form two radially dis- 
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tinct population of RR Lyrae stars through a combina- 
tion of hierarchical mergers which could have been dis- 
sipative and dissipat ionless (e.g., iBekki &: Chibal 120011 : 
IChiba &; Beers! [2000t ). It is possible that systems ac- 
creted at early times contained considerable gas, and 
would have resulted in dissipative merging. However, 
systems accreted at later times would have contained less 
gas due to star formation and ejection from supernovae 
explosions. This late-time merging would have resulted 
in dissipatio nless merging . Using a large sample of stel- 
lar spectra, iCarollo et alj (|2007f ) found that the stellar 
halo is composed of at least two spatially, kinematically, 
and chemica lly distinct compone nts, which agrees with 
our results. ICarollo et al.l (|2007l ) argue that their data 
can be explained solely within the context of hierarchi- 
cal merging. This scenario could account for the two 
spatially distinct components in the LONEOS - 1 sample 
as well. 

We have found that the stellar halo is predominately 
Ool, with a significant OoII component. We cannot 
rule out whether the stellar halo c ontai ns an interme- 
diate Oosterhoff class. ISieeell (|2006l ) and lDall'Ora et all 
(2006) have found that the recently discovered Bootes 
dSph galaxy to be of OoII type. Since the Bootes dSph 
is metal-poor and is an OoII, it supports the argument 
that the Oosterhoff effect is a continuous effect which 
is predominately associated with metallicity. This new 
finding raises the possibility the OoII component of the 
Galactic stellar halo was formed from accreted systems 
like Bootes. 

The chemical abundance properties of the stellar halo 
and the dwarf galaxies are distinct (e.g., IVenn et abl 
2004). The stellar halo is more metal-poor, but with en- 
hanced \ a/Fe] abundance s , than the Local Group dwarf 
galaxies. iRobertson et al.l (|2005| ) have shown that such a 
discrepancy can be explained within the CDM paradigm. 
Using cosmologically-motivated star formation and gas 
accretion histories, they argue that most of the stellar 
halo was formed from short-lived satellites which were 
accreted and destroyed at early times. Since these early 
progenitor satellites were short-lived, there is limited 
time for stellar enrichment and the dominant contrib- 
utor to the enrichment is from Type II supernovae (i.e., 
a-enrichment) . However, the present day Local Group 
dwarf galaxies have had a long time to produce metals 
from both Type I and II supernovae. The bottom line 
is that the dwarf galaxies that we observe in the Local 
Group may not be representative of the progenitors of 
the stellar halo. A detailed study of metallicity of the 
two Oosterhoff halo components will allow us determine 
if there is an age difference between these two Oo com- 
ponents. 

6. CONCLUSIONS 

We present a catalog of 838 aft-type RR Lyrae in the 
Galactic stellar halo from the LONEOS - I survey. We 
find evidence for two distinct components in this sample. 
We associate these two components with a manifestation 
of the Oosterhoff effect in the stellar halo. In addition, we 
find that these two components have significantly differ- 
ent radial distributions. Without flattening, the Ool has 
a power-law exponent of —2.26 ±0.1 while the OoII com- 
ponent has a significantly steeper slope with — 2.88±0.04. 
The flattening of the halo cannot be constrained by this 



data set since the distribution on the sky is limited. We 
performed numerous tests and simulations to verify the 
existence of the two components and their difference in 
their Galactocentric radial distributions. Our finding is 
robust with respect to potential systematic errors due to 
e.g. detection efficiencies, zeropoint problems, color dif- 
ferences, camera systematics, etc. It is likely that the 
OoII component potentially represents an older popu- 
lation and traces an early dissipative collapse or early 
dissipative merging activity. The Ool component repre- 
sents the later accretion of younger systems. Spectro- 
scopic data would be very valuable, the radial velocities 
will tell us if these Ool and OoII components have dif- 
ferent kinematic properties and occupy different regions 
in phase-space. Spectroscopic data could also supply ev- 
idence for chemical differences. Chemical differences in 
these two populations will give us information regarding 
the environment and age of these populations. 

In 2000, the LONEOS camera (LONEOS - II ) was 
upgraded, and now consists of two 2k x 4k backside- 
illuminated CCDs, giving it an 8.1 square degree FOV 
(2.8 arcsec/pixel) , and reaches to a depth of R ~ 19.3. 
Approximately one quarter of the sky is covered with 
more than 40 epochs of measurements. Photometric re- 
duction and analysis of the new camera data has com- 
menced, and results will be discussed in future papers. 
The LONEOS - II data should also increase our RR 
Lyrae sample by an order of magnitude with a wider 
and deeper survey area. The data from the LONEOS - 
II camera will produce cleaner period-amplitude distri- 
butions. This will be enable us to search for additional 
"fine structure" in the period-amplitude distribution. 
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